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Abstract 



In this letter we investigate the deep Newtonian regime of the 
MOND paradigm from a purely phenomenological point of view by 
exploiting the least-square estimated corrections to the secular rates 
of the perihelia of the inner and of some of the outer planets of the 
Solar System by E.V. Pitjeva with the EPM2004 ephemerides. By 
using fj.(x) w 1 — fco(l/a;)" for the interpolating MONDian function, 
and by assuming that feg , considered body-independent so to avoid vi- 
olations of the equivalence principle, experiences no spatial variations 
throughout the Solar System we tightly constrain n with the ratios of 
the perihelion precessions for different pairs of planets. We find that 
the range 1 < n < 2 is neatly excluded at much more than 3 — tr level. 
Such a test would greatly benefit from the use of extra-precessions of 
perihelia independently estimated by other groups as well. 
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1 Introduction 

The MOND scheme (MOdified Newtonian Dynamics) was put forth by Mil- 
grom in Ref. [l] in order to phenomenologically explain two basic facts con- 
cerning spiral galaxies without resorting to the concept of hidden dark mat- 
ter: the asymptotic flatness of the rotation curves of spiral galaxies p] and 
the Tully-Fisher law which is a well-defined relationship between the rota- 
tion velocity in spiral galaxies and their luminosity [3j. 

Viewed as a modification of gravitj0, MOND predicts that the gravita- 
tional acceleration Ag felt by a particle in the field of a distribution of mass 



^It can also be considered as a modification of tlie inertia of a particle under the action 
of a generic force F. 
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where An is the Newtonian acceleration, Aq is an acceleration scale which 
different, independent ensembles of observations set to [3] Aq = 1.2 x lO^^'^ 
m s~^, and //(x) is an interpolating function which approximates 1 for x ^ 1, 
i.e. for accelerations larger than Aq] for x ^ 1 = x, so that in such a 
strongly MONDian regime Ag ~ A-^Aq. 

Here we wish to investigate the deep Newtonian regime (x 2> 1) in 
view of recent advances in planetary orbit determination occurred for the 
inner planets of the Solar System. For a quite general class of interpolating 
functions, /i(x) can be cast into the form [T] 

//(x) « 1 - /cq ( i 

which yields a modified gravitational acceleration [5] 



Ag PS An 



(3) 



Note that the most commonly used expressions for /x(x), i.e. [T] 



/^(a;) = J-^ ^ (4) 
Vl + x^ 



and [6] 



= YT-x^ (5) 

reduce to eq. ([2]) for = 1/2, n = 2 and fco = l,n = 1, respectively, in the 
appropriate limit. For a recent review of many aspects of MOND as various 
attempts to theoretically justify it see Ref. [7] and Ref. 



2 Constraints from planetary perihelion preces- 
sions 

It can be shown that eq. ([3]) affects the orbital motion of a test particle 
in the field of a central mass M with a secular rate of the longitude of the 
pericenter [U] 

w = ^ a^" 2 + O(e^) Qa^ (6) 
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where ru = \/ GM/Aq, a and e are the orbit's semimajor axis and eccen- 
tricity, respectively, and 



konVGM ^ 3 
Q= ' z = 2n--. 



(7) 



The expression of eq. ([6]) can be fruitfuhy used in conjunction with the 
corrections to the known Netonian/Einsteinian secular rates of the perihe- 
lia of the inner planets of the Solar System phenomenologically estimated 
as least-square solve-for parameters in Ref. [TO] by fitting a planetary data 
set spanning almost one century with the dynamical force models of the 
EPM2004 ephemerides [llj . The same procedure was followed also for 
some of the outer planets [12]. Since such models fully include Newtonian 
and Einsteinian gravity, such estimated extra-precessions account, in prin- 
ciple, for all unmodelled physical effects possibly present in nature like, e.g., 
MOND. If and when other groups will independently estimate their own 
extra-precessions of perihelia it will be possible to enlarge and enforce the 
present test. For a search of other MOND-like effects in the Solar System 
see Ref. [13], while the possibility of testing a MONDian violation of the 
Newton's second law in a terrestrial environment is discussed in Ref. [14j . 
By assuming that MOND does not violate the equivalence principle, i.e. 
and n are not body-dependent, it is possible to consider for a generic pair 
of planets A and B the ratio of their perihelion rates getting 
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By defining 



and 



n 



l(B) 



■w 



w 



(B) 



j(A) 
7P) 



it is possible to construct 



n-e. 



(8) 

(9) 
(10) 
(11) 



"^It is implicitly assumed that fco does not experience spatial variations, according to the 
MONDian point of view for which modifications of Newtonian gravity does not depend 
on distance but on acceleration only 4 . Of course, our test based on the ratio of the 
perihelia is vaUd for the case fco 7^ 0, n 7^ 0. 
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if, for a given value of n, the quantity |r„|, computed with the extra-rates 
of periheha of Ref. [10] estimated without including any exotic acceleration 
with respect to standard Newton-Einstein one in the suite of the dynamical 
force models used to fit the data, turns out to be incompatible with zero 
within the errors, i.e. if |r„|/(5r„ > 1, that value of n must be discarded. 
Note that our test makes sense for ko ^ 0, as it is just the case from galactic 
data. The uncertainty in r„ can be conservatively assessed as 



5r„ <6U + 5G„, (12) 

with 

5n < ini 



|zi7(A)| |tZ7(B)| 



6Qn < zQn 



(13) 



(14) 



The linear sum of the individual errors in eq. (jl3p accounts for the exist- 
ing correlations among the estimated perihelia corrections, which reach a 
maximum of about 20% for Mercury and the Earth (Pitjeva, private com- 
munication 2005). 

By choosing A=Mars and B=Mercury, from Table [1] we get Figured] in 
which we plot |r|/(5r for 1 < n < 2. It turns out that the corrections of 
order 0{e^) to eq. ([6]) are negligible in the sense that their inclusion in the 
calculation does not alter the results for |r|/(5r. Moreover, 6Q is far smaller 
than (5n even by re-scaling the formal errors in the semimajor axes by a 
factor 10 or more. As can be noted, |r| is always incompatible with zero at 
much more than 3 — cr level, thus ruling out the interval [1,2] for n. Figure 
[2] refers to A=Earth, B=Mercury: it yields the same conclusions. It can 
be shown that the same holds also for A=Mars, B=Venus and A=Earth, 
B=Venus. It maybe interesting to note that, although the errors in the 
perihelion rates quoted in Table [1] are not the mere, formal ones, should 
one decide to re-scale them by a factor 10 our conclusions would remain 
unchanged, apart from the number of a which would pass from about 180 
to 18 (A=Mars, B=Mercury), or from 60 to 6 (A=Earth, B=Mercury) for 
n = 2. 

Another way to tackle the problem is to construct a x^-like quantity ^ 
defined as follows. The role of the observables Ok is played by the ratios of 
the perihelia precessions 11^. for all the k = 1,2, ...30 pairs A/B and B/A of 
planets, including also the gaseous giant ones whose data are quoted in Table 
[T] The computed, or predicted, quantities Ck are the ratios (aA/aB)^""^ for 
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Figure 1: |r|/(5r for Mars and Mercury and 1 < n < 2. As can be noted, 
|r| 7^ for all values of n. 
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Figure 2: |r|/(5r for the Earth and Mercury and 1 < n < 2. As can be 
noted, |r| 7^ for ah values of n. 
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the corresponding pairs of planets. The errors cr^ are 511^ because of the 
neghgible impact of the uncertainties of the semimajor axes, while we will 
assume the number of different pairs constructed, i.e. 30, for the number of 
degrees of freedom d. Thus, we obtain 



which largely confirms our previous conclusion. It turns out that ^ ^ 1 also 
for 1 < n < 2. 

3 Discussion and conclusions 

A major outcome of our analysis is that n = 2 is neatly ruled out indepen- 
dently of ko and contrary to what obtained in Ref. [9j in which only 
the perihelion of Mars was used by keeping fixed /cq ~ 1 (and using the 
commonly accepted value [4] = 1.2 x lO"^'^ m s~^). The results of the 
analysis presented here seem to suggest that a reconsideration of the match- 
ing between the deep Newtonian and MONDian regimes should be looked 
for. However, caution is in order because the present analysis is based upon 
the extra-precessions of perihelia estimated by only one team; it would be 
of great importance to have at disposal corrections to the perihelia rates 
determined independently by other groups as well. 
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Table 1: Semimajor axes a, in AU (1 AU= 1.49597870691 x 10^^ m), and 
phenomenologically estimated corrections to the Newton-Einstein perihelion 
secular rates, in arcseconds per century (" cy^^), of the inner [10] and some 
of the outer [12] planets. Also the associated errors are quoted: they are 
in m for a (see Ref. [U]) and in " cy~^ for w [lOj [12]. For the semimajor 
axes they are the formal, statistical ones, while for the perihelia of the 
inner planets they are realistic in the sense that they were obtained from 
comparison of many different solutions with different sets of parameters and 
observations (Pitjeva, private communication 2005). The errors quoted here 
for the perihelia of the outer planets are the formal ones re-scaled by a factor 
10 in order to get realistic estimates for them. 



Planet 


a (AU) 


5a (m) 


w {" cy-^) 


i" cy-i) 


Mercury 


0.38709893 


0.105 


-0.0036 


0.0050 


Earth 


1.00000011 


0.146 


-0.0002 


0.0004 


Mars 


1.52366231 


0.657 


0.0001 


0.0005 


Jupiter 


5.20336301 


639 


0.0062 


0.036 


Saturn 


9.53707032 


4,222 


-0.92 


2.9 


Uranus 


19.19126393 


38,484 


0.57 
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